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Abstract 

A mean field model is presented for the configuration dependent effective demagnetizing and 
anisotropy fields in assemblies of exchange decoupled magnetic particles of arbitrary shape which 
are expressed in terms of the demagnetizing factors of the particles and the volumetric shape con- 
taining the assembly. Perpendicularly magnetized 2D assemblies have been considered, for which 
it is shown that the demagnetizing field is lower than the continuous thin film. As an example 
of these 2D systems, arrays of bistable cylindrical nanowires have been characterized by rema- 
nence curves as well as ferromagnetic resonance, which have served to show the correspondence 
of these measurements with the model and also to validate the mean field approach. Linear chains 
of cylinders and spheres have been analyzed leading to simple expressions to describe the easy 
axis rotation induced by the interaction field in chains of low aspect ratio cylindrical particles, and 
the dipolar magnetic anisotropy observed in the linear chain of spheres. These examples serve to 
underline the dependence on the dipolar interaction field and effective demagnetizing factor on the 
contributions that arise from the shape of the outer volume. 
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Assemblies or arrays of exchange decoupled magnets represent a very large class of magnetic 
systems which are of great interest in different subjects and problems in magnetism at both, fun- 
damental and technological level." These include assemblies in which the particles are arranged 
or dispersed on a substrate or embedded in different media such as liquids, polymers or other non 
magnetic materials.-- Depending on their spatial disposition, and thus the shape or outer volume 
that contains the assembly, they form one, two or 3D heterostructures or composites 

The magnetic properties of a single isolated particle are mainly governed by the shape and 
magnetocrystalline anisotropy and eventually a magnetoelastic contribution. In an assembly of 
many of such particles, their individual properties are further modified by the interaction field 
produced by the other particles. If the particles are not allowed to come into contact, then the 
interaction is purely dipolar and has the same magnetostatic origin as the demagnetizing field 
and shape anisotropy associated to each particle. In this sense, from a mean field perspective, 
an effective demagnetizing field which includes the particle self demagnetizing field as well as 
the dipolar interaction is a common feature of any assembly of exchange decoupled assembly 
of particles.-ii^ Consequently, these considerations should apply to any assembly of particles, 
whether the particles are magnetically hard or soft, and regardless of the approach or level of 
approximation used to describe them. In particular, it should be possible to interpret the vast 
phenomenology of effects associated to the dipolar interaction or the effective demagnetizing field, 
within this framework.— — However, this is usually complicated by the difficulty of finding either 
adequate expressions or reliable measurements for the dipolar interaction. Furthermore, extensions 
to include the magnetization or configuration dependency of the dipolar interaction field adds 
additional difficulties for establishing a suitable description for these systems. 

A coherent formulation is required for the effective demagnetizing field which allows to include 
the configuration dependent interaction field for an assembly of particles of arbitrary shape and 
outer volume and that allows determining how the interaction field will modify the demagnetizing 
field, the effective anisotropy or total energy and coercivity of the system. 

The present study is focused in deriving mean field expressions for the configuration depen- 
dent effective demagnetizing and anisotropy fields in assemblies of exchange decoupled magnetic 
particles, which depend on both the shape of the particles, the volume containing them and the 
packing fraction. The magnetostatic properties of 2D and ID assemblies have been analyzed and 
the results show that the dipolar interaction depends not only on the interparticle distance, but 
also on the outer shape of the assembly. In particular it is shown that this outer shape contri- 
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bution to the dipolar field plays an important role in the coercivity of the assembly, the limiting 
values of the effective demagnetizing fields and the magnetization reorientation transition due to 
the dipolar interaction. Remanence curves and ferromagnetic resonance measurements done on 
arrays of bistable cylindrical nanowires have served first to identify the correspondence of these 
measurements with the model, and second, to validate this mean field approach. 

I. EXPERIMENTAL 

Arrays of Co, NigiFcig and Co55Fe45 nanowires having a low density (low porosity P) and 
small diameter (29 > > 40 nm) have been grown by electrodeposition into the pores of 21/im 
thick lab-made track-etched polycarbonate (PC) membranes, in which the pores are parallel to 
each other but randomly distributed.^" Also, two samples have been fabricated using anodized 
allumina templates in order to have higher porosities (10 and 15%). Full details of the preparation 
process can be found elsewhere.— 

For the electrodeposition a Cr/Au layer is evaporated previously on one side of the membrane 
to serve as a cathode and deposition is done at a constant potential using a Ag/AgCl reference 
electrode. For CoFe a 40 g/1 FeS04 -i- 80 g/1 C0SO4 -1- 30 g/1 H3BO3 electrolyte was used with a 
potential of V=-0.9 V, while for NiFe the electrolyte contained 5.6 g/1 FeS04 -1- 131.4 g/1 NiS04 -1- 
30 g/1 H3BO3 and deposition is done at V"=-l.l V. Cobalt nanowires have been grown at V=-IY 
using a 238.5 g/1 C0SO4 -1- 30 g/1 H3BO3 electrolyte with the pH set to 2.0 by addition of H2SO4 to 
favor a polycrystalline fcc-like Co structure with no magnetocrystalline anisotropy contribution.— 
Table IJ shows the details of the samples considered. 

Ferromagnetic resonance (FMR) measurements have been done in the field swept mode with 
frequencies ranging from 1 up to 50 GHz using a 150 /im wide micro stripline with the DC 
magnetic field applied parallel to the long axis of the wires, as detailed elsewhere.— "^^^ For all the 
samples considered, the transmission spectra is recorded at a fixed frequency while the magnetic 
field is swept from 10 down to kOe. These measurements are repeated for different frequencies 
and from the collection of frequencies and their respective resonance field, the dispersion relation 
is obtained. Figure[I](a) shows typical dispersion relations obtained on the samples. The dispersion 
relation for an array of infinitely long cylindrical nanowires with the field applied parallel to the 
wires which are in the saturated state 

f = j{HR,, + H,fj), (1) 
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TABLE I. Co, NiFe and CoFe nanowire samples used in this study numbered as si to si 6. For each sample 
the following quantities are given: the wire diameter <I>, nominal packing fraction P/v> the FMR effective 
field i^e//» the dipolar interaction coefficient a^, and the determined values of the packing fraction Pm and 
saturation magnetization M*. 
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where Hr^s and H^ff are the resonance and effective anisotropy field, respectively. From the 
linear fit of the dispersion relation, see Fig. [11(a); the values of H^ff have been determined for 
each sample, and the corresponding values are given in Table ID 

Isothermal Remanence Magnetization (IRM) and DC demagnetization (DCD) remanence 
curves measurements with the field applied parallel to the wires axis, have also been done in 
all the samples, from which the interaction field coefficient along the wire axis, a^, has been 
determined as2^ 

a, = 2(i70-^-i70), (2) 
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FIG. 1. (a) Typical dispersion relations obtained on anays of NWs with the field applied parallel to the 
wires, the measurements ai^e shown as symbols, while the continuous line is the best fit using Eq. ([T|). (b) 
IRM (blue) and DCD (red) remanence curves measured on sample SIO [NiFe] with the field applied parallel 
to the NWs axis. 



where is the field value at which the normalized IRM curve is 0.5, while is the field value 
at which the DCD curve is zero, as shown in figure[l](b). The values determined for the interaction 
field coefficient are given in table H 



n. THE EFFECTIVE DEMAGNETIZING FIELD 

Consider an assembly of exchange decoupled identical particles where each one has a de- 
magnetizing factor given by N, called the inner demagnetizing factor, contained within a macro- 
scopic volume described by an outer demagnetizing factor , as schematically shown in fig- 
ure |2l Furthermore, all the particles are aligned so their easy axis point in the same direction. 
The energy density of a given particle can include the magnetocrystalline {Kmc), shape {Ks) 
and magneto elastic (Kme) anisotropics, which determine the total anisotropy of the particle, 
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= Kmc + Ks + Kme, as well as a contribution due to the interparticle interaction {Eint) 

E = KAsin^e + Eint (3) 



Where the shape anisotropy and the dipolar interaction are of magnetostatic origin and are related 
to the effective demagnetizing fields, which are considered in the following sections. Furthermore, 
in the following it will be assumed that the particles are homogeneously magnetized regardless of 
their shape. 

A. Saturated state 

The effective or total demagnetizing factor can be expressed using the interpolation procedure 
introduced by Netzelmann.— To simplify the notation, in the following, the effective or total de- 
magnetizing factor (Neff) will be denoted as Nt which can be expressed in terms of the packing 
fraction P asr^^ 

Nt = (l- P)N + PN^. (4) 

However, for the treatment of an assembly of magnetic particles it is more useful to rewrite this 
last expression as, 

Nt = N + {N+ - N)P. (5) 

The first term on the right side corresponds to the self demagnetizing factor of the particles that 
makeup the assembly, this is, Ns^if, while the second term contains all the contributions that 
depend on the packing and corresponds to the dipolar interaction contribution, Ndip, then 

Nt = Nseif + Ndip. (6) 

When P = the demagnetizing factor reduces to that of the isolated non-interacting particle. 
While on the opposite limit, P — )■ 1 this leads to the outer demagnetizing factor. If both the inner 
and outer demagnetizing factors are diagonal and their trace is Tr(A^) = An, then from Eq. (|5]l, 

Tt{Nt) = Tr(iV) + [Tr(iV+) - Tr(iV)]P, (7) 

and from the form of the following properties of Nt are obtained, 

Tr(iVT) = 47r, (8) 
Tr{N,,p) = 0. (9) 
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FIG. 2. Schematics of an assembly of identical particles characterized by an inner demagnetizing factor N 
contained within a macroscopic volume geometrically characterized by the outer demagnetizing factor 
(dashed lines). 

The z-th component of the effective demagnetizing field in the saturated state for an assembly of 
particles follows from equation (|5]), 

H^' = M,N, + (iV+ - N,)M,P, (10) 

where the first term is the usual self demagnetizing field of the particle, Hf = MgNi, and the 
second term is the dipolar interaction field, this is, 

Hi^ = {N+-N,)M,P, (11) 

The total or effective magnetostatic anisotropy field in the saturated state defined as = 
Ms{N^ — N^), where the hard axis is along the x-direction and the easy axis is taken along 
the 2-axis, then, 

= MsAN + (AiV+ - AN)MsP, (12) 

The first term on the right side is the shape anisotropy field of the individual particles (Hs), while 
the second term corresponds to the total or effective dipolar field in the saturated state, 

Hdip = {AN+ -AN)MsP, (13) 

The total shape (magnetostatic) anisotropy energy is given by = MgHj /2, which includes 
both the shape anisotropy of the particle and the contribution of the dipolar interaction. This 
anisotropy constant now regroups Ks and Eint in equation ([3]). 

Finally, from the expressions obtained for the interaction field there are two other relations 
of interest that follow when both the inner and outer volumes have in plane symmetry, so that 
Nx = Ny. Due to the symmetry, the in-plane components of the interaction field are equal, 
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Hdip = -^dip' ^'^d since Tr{Ndip)=0, according to equation the following relation between the 
components of the dipolar fields is obtained, 

Hl,^ = -2Hl^. (14) 

This shows that when both the inner and outer volume have rotational symmetry, the dipolar 
interaction field along the symmetry axis is twice the dipolar field in the hard axis with opposite 
sign. This, for example, has been obtained by several authors .■^^ A consequence of this result is 
that the dipolar part of the total or effective anisotropy field {Hdip = H^^^ — H^-^) is 
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Hdip = --H^ip, (15) 



which provides a relation between the component of the interaction field along the easy axis, with 
the net contribution of the interaction field to the total anisotropy field. 

B. Non saturated states 

To extend the previous expressions to the non-saturated case, consider that the particles are 
bistable, this is, they can only be magnetized along their easy axis in both the positive or negative 
direction. The normalized magnetization, m = M{H) / Ms, so that —1 < m < 1. Any magnetic 
state called hereafter m, can be written in terms of the fraction of particles magnetized in the 
positive (m+) and negative (m„) directions, with < m± < 1, and m = m+ — m„. Moreover, 
since the number of particles is constant, m+ + m_ = 1, and m± can be obtained from the value 
of m as, 

m± = . (16) 

A recent FMR study done on non saturated arrays of bistable NWs lead to the expressions for the 
effective field, the dipolar interaction term and the FMR dispersion relation, that include explicitly 
the magnetic configuration of the system. ^^-^^ In particular, it was shown that to extend the known 
expressions for the saturated case to the non-saturated case, it is necessary to rewrite the dipolar 
interaction term as a parametric expression of m±r^^ Assuming a dipolar interaction field of the 
form am, it implies that the dependence of the interaction field on m is not as a direct product but 
instead, requires using Eq. (fT6] ). this is, 

Hdipm± = ± ^^m, (17) 



which at saturation reduces to the expected value. The ± sign is introduced to include both positive 
(m = 1) and negative (m = — 1) saturation. 

Using Eq. (fTOl ). the ith-component of the effective demagnetizing field is, 

H^' = MsN, + {Nt -N,)^± (iV+ - N,)^m. (18) 

The dipolar contribution now contains two terms, the first one is constant, while the second one is 
proportional to m and corresponds to the magnetization dependent part of the demagnetizing field. 
If the magnetization dependent interaction field is taken as am, then from the last term on the right 
hand side, the i-th component of the dipolar interaction field coefficient a can be identified as, 

a^ = ^ = {Nt-m^. (19) 
Then equation (fT8l) . can be written as, 

= H^ + a,± aim, (20) 

where Hf^ is the demagnetizing field of the isolated particle. From Eq. (fT2l) . the effective or total 
magneto static anisotropy field is, 

= Hs + ar ± arm, (21) 



and ax = a^ — a^ is the total dipolar field coefficient. At saturation m = 1, both equations ([20] 
and (|2T]) reduce to Eqs. (fTOl l and (fT2)) . respectively. 



III. APPLICATIONS 



A. Two dimensional arrays with rotational symmetry 

An important class of magnetic assemblies are two dimensional arrays of exchange decoupled 
nanoparticles, ideally a monolayer of single domain particles, with perpendicular magnetization. 
These correspond, for example, to systems obtained by lithography, including perpendicular bit 
pattemed media, also granular thin films with columnar structure as those used for perpendicular 
recording media, self assembled monolayers and nanowire arrays.-^ 

For these systems one can assume that the height of the particles is very small compared with 
the lateral dimensions of the entire array and the volume containing the particles can be considered 
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as an infinite thin film, so that = Ny = and A^^ = An, as depicted in figure |2] (a) and (b). 
Using these values in Eq. (flOl) and assuming particles with in-plane symmetry so that = Ny 
and using = 4:1t — 2Nx to express quantities in terms of N^, the effective demagnetizing field 
perpendicular to the plane {z) is 

= A-nMs -2N^{1- P)Ms. (22) 

Which shows that the effective demagnetizing field of a thin film made of exchange decoupled 
entities is less than Air Ms by a quantity equal to 2Nx{l — P)Ms. 
The total shape anisotropy field [Eq. (|2TI )1 is, 

= (3iV, - 47r)M, - ^N^M^P T ^N^MsPm, (23) 

the first term on the right side is the shape anisotropy of the particle, while the second and third 
term correspond to the dipolar part of the total anisotropy field, for m = 1, 

Hmp = -^N^MsP, (24) 

which is a particular case of Eq. (fTST ). As seen from this expression H^ip is only a function of 
and P. Since the particles have in-plane symmetry, then = An — 2N^^, from where it follows 
that if < < An, then < A^^ < 2n. And from the previous expression for Hdip, it is 
possible to establish an upper bound for the dipolar part of the total or effective anisotropy. Taking 

N, = 2n, 

H^;^ = -QnMsP (25) 

The dipolar contribution to the total anisotropy field is antiferromagnetic, inferred by the negative 
sign in Eq. (|24|) . as expected for a 2D array of nanomagnets with perpendicular anisotropy. 

B. Cylindrical Nanowires 

In the particular case of a 2D array of circular cylinders of arbitrary height aligned parallel to 
each other so their axes are along the z axis, as shown in figure|2](b), their effective demagnetizing 
fields along the easy and hard directions are given by equation (flOl ). As in the previous section, 
the outer demagnetizing factor is taken as an infinite thin film. For a cylinder of arbitrary aspect 
ratio (wire height divided by its diameter), can be determined numerically.^ However, for the 
NWs considered in this study, the wire height is typically of the order of 20/im so the aspect ratio 
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is very high, and they can be considered as infinite so their demagnetizing factors are A^^ = and 
Nx = Ny = 2n and from equation (fTOl ) the demagnetizing field along the easy and hard axis are, 

H^' = AttMsP, (26) 
H^'' = 2'kMs - 2ttMsP. (27) 

Due to the symmetry properties of both inner and outer demagnetizing factors, H;^^ = h!^^ and 
the dipolar terms satisfy both equations ^ and (fT4] ). The effective anisotropy field is, 

= 2-kMs - GtiM.P, (28) 

that corresponds to the known expression for the effective field for an array of infinite cylindrical 
nanowires in the saturated state.— The effective dipolar interaction field is HfHp = —GttMsP, and 
using Eqs. (fT9l ) and (|2T]) . the magnetization dependent part of the effective interaction field is, 

aT = = 3nM,P, (29) 



in agreement with the expression obtained in Ref. [|29ll . 

On the other hand, the measurement of the interaction field using Eq. Q and the IRM and 
DCD remanence curves are done solely along the wire axis, this is, the measurement is related to 



the component of the interaction field along this direction, as pointed out in Ref. 123|. From Eq. 
(|26|) . and Eqs. (fT9l) and (l20l) . the component of the magnetization dependent part of the interaction 
field along the wire axis is 

= 2nMsP. (30) 

Then, Eqs. (l28l) . (|29|) and (l30l) provide different expressions for the interaction field which have 
different physical meaning. Moreover, depending on the measuring technique, each of these three 
quantities can be measured independently ,■^^'2^ so it is important to distinguish between them. 

Since the effective anisotropy field, can be determined from the dispersion relation obtained 
from the FMR measurements, using equation ©, then from Eq. (l28l) . HfHp = H^ — 2'kMs. So H^ip 
has been determined for each sample from the FMR measurements using the nominal values of Mg 
for Co (1400 emu/cm^), NigiFcig (800 emu/cm^), and Co55Fe45 (1900 emu/cm^) and the values 
of Heff given in Table IB These values where then divided by 3 so they correspond numerically to 
Eq. (l30l) which is the value of the component of the interaction field along the wire axis measured 
with the IRM and DCD remanence curves. Figure [3] shows the values obtained using the IRM and 
DCD remanence curves, plotted as a function of the values obtained by FMR, Hdip/S. 
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FIG. 3. Easy axis component of the interaction field determined by magnetometry, = iTrMgP, plotted as 
a function of Hdip/S, where Hdip = GirMsP is the total interaction field measured by FMR in the saturated 
state, for all the nanowire samples listed in Table U 



These results show a very good agreement, which shows that both Eqs. (1281) and (1301) provide 
correct results for the interaction field. Moreover, these two expressions provide expressions for 
quantities that can be measured independently, so its possible to find a self consistent method 
to determine both Ms and P. Since Hj^ and az are determined from FMR and magnetometry 
measurements, respectively, then combining Eqs. (|28] ) and (|30l) leads to the following expressions 
for M, and P. 



+ Mz 

Ms = (31) 

- (32) 



The values of and measured by FMR and the remanence curves, respectively, that are 
given in Table H] have been used as input in Eqs. (|3T| ) and (|32l ) to obtain the corresponding values 
of Ms and P. The results are given in Table |I] in the columns labeled Pm. and M*, respectively. 
As a first point, the values of M* show a very good agreement with the known values for Co, 
NiFe and CoFe. Furthermore the values for NiFe and CoFe alloys are in good agreement with 
those determined solely by FMR and reported elsewhere.-22i^ Regarding the values of the template 
porosity, comparing the nominal values determined by scanning electron microscopy P^r with 
those determined by the measurements P^, a very good agreement is also found in practically all 
the samples. This procedure is general and can be extended to other assemblies by solving for Mg 
and P using Eqs. ([H]) and (fT9l) . 
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C. Linear chain 



So far, only 2D systems with perpendicular magnetization have been considered. Another 
interesting example is that of the linear chain of particles and particularly for circular cylinders 
of different aspect ratio and spheres, as those shown in figure [2] (c) and (d), respectively. For an 
infinite number of particles, the outer demagnetizing factor is an infinite cylinder of circular cross 
section whose axis is assumed along the z axis so A'^^ = = 2tx and = 0. While for 
both spheres and cylinders, the components of the inner demagnetizing factor are = Ny and 

= 4:TV — 2Nx. From equation (fT2] ). the effective anisotropy field is, 

= M, (iV, - N,) + M, (67r - 3iV,)P. (33) 

For the calculations, the demagnetizing factors and the packing fraction are required. For the chain 
of spheres, = Ny = = Att/3 and the packing fraction is P = {2/3)(f)/d where d is the 
center to center distance and (p is the diameter of the sphere. For the cylinders, is determined 
using known expressions as a function of the aspect ratio r = h/cp where h is the height of 
the cylinder,— and in this case, P = h/d. Figured (a) shows the reduced effective anisotropy 
AA^ = H^/Ms for linear chains of cylinders with different aspect ratio (r=0.5, 0.8, 2.5 and 5) as 
well as for a chain of spheres (dashed line). 

Consider first the case of the linear chain of cylinders. For a single, isolated and non interacting 
circular cylinder, the critical aspect ratio is r = 0.9065 above this value, N,j. > and the easy 
axis is along the cylinder axis and at lower values < and it is perpendicular to the axis. As 
seen in Fig. |4](a) at large distances the anisotropy tends asymptotically to the expected value of 
the isolated non-interacting cylinder, which yield negative values for r =0.5 and 0.8, and positive 
for T =2.5 and 5. 

As the particles are brought closer the interaction increases. For r > 0.9065, the interaction is 
ferromagnetic as it favors head to tail alignment of the magnetization, which results in an increase 
of the effective anisotropy as the packing fraction increases. For low aspect ratios (r < 0.9065) 
the easy axis of a given cylinder is perpendicular to the chain axis and the first term in Eq. (l33T l 
is negative, so the interaction becomes antiferromagnetic as it has the opposite sign, and com- 
petes against it. As the packing fraction is increased, the value of the interaction increases and 
overcomes the anisotropy of the particle leading to a reorientation of the magnetization easy axis. 

In other words, to have a reorientation of the magnetization easy axis due to the dipolar inter- 
action, the interaction has to be of the demagnetizing type (antiferromagnetic) with respect to the 
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FIG. 4. (a) Reduced effective anisotropy as a function of the reduced center to center distance d/cj) fox linear 
chains of cylinders with aspect ratio r=0.5, 0.8, 2.5 and 5 as well as for spheres (dashed line), (b) Zero 
anisotropy curve for a linear chain of cylindrical particles of aspect ratio r as a function of the reduced 
center to center distance. 

easy axis direction of the particle when considered isolated. For the linear chain of cylinders there 
is a critical aspect ratio for the easy axis reversal which depends on the inter particle distance, that 
correspond to those values for which the effective anisotropy vanishes. Setting = in Eq. 
(1331 ) and using 2Nx = An — N^, the critical packing fraction can be expressed as, 

which is only valid for A^^ > 47r/3. Here A^^ = Nz{t) gives the dependence on the aspect ratio 
of the cylinder which is given by well known expressions.— Figure |4](b) shows the critical aspect 
ratio variation with the reduced center to center distance for a linear chain of cylinders. Above 
this curve, the magnetization lies along the cylinder axis and the interaction between them is 
ferromagnetic. Below this line, the magnetization easy axis is perpendicular to the cylinder axis 
and the dipolar interaction is antiferromagnetic. At large distances, this curve tends asymptotically 
to r = 0.9065 (horizontal dashed line) which is the critical aspect ratio for an isolated non- 
interacting cylinder with A'^ = 47r/3.— The highest packing fraction required to reverse the easy 
axis is obtained when the axial demagnetizing factor of the disk tends to 47r, which from Eq. 
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(I34l) . corresponds to Pc = 2/3. 

Finally, the linear chain of spheres presents features that are interesting as seen in Fig. |4](a), in 
particular the existence of an effective anisotropy and its increase as the packing fraction increases, 
which in light of Eq. (l33l) is entirely due to the dipolar interaction. Indeed, as seen in the figure, 
the effective anisotropy tends to zero as the distance between particles is increased consistent with 
a single sphere with no shape anisotropy. Moreover, the maximum value of the anisotropy is 
attained when the spheres come in contact and this value is lower than 27r, which shows that there 
is not a full equivalence between a cylindrical wire and the chain of spheres. This point has been 
brought to light recently and its important since modeling of cylindrical wires is commonly done 
by considering them as a chain of spheres.-^^ From Eq. (l33l) . the maximum value of the effective 
anisotropy for the chain of spheres is reached at P = 2/3 and, 

HLa. = Y^s, (35) 

as indicated in Fig. |4](a). 
IV. DISCUSSION 

Expressions for the magnetization dependent demagnetizing field and the magnetostatic anisotropy 
field (energy) have been derived. These are general as the inner and outer volume are arbitrary and 
apply to assemblies of identical exchange-decoupled, bistable particles with a common easy axis. 
Moreover, this approach is a mean field approximation which only accounts for magnetostatic ef- 
fects. In consequence, all the particles are equivalent which is only valid if the assembly contains 
a very large number of particles. 

Regarding the form of the effective demagnetizing field in the saturated state, Eq. (flOl) . it 
contains the self demagnetizing field of the individual particle and the component of the dipo- 
lar interaction field, Eq. (fTTT ) which is consistent with previous expressions obtained for nearly 
spherical particles.-'* On the other hand, when compared to other expressions,— which provide 
non physical results for particular limiting cases,— Eq. (flOl ) provides correct limiting values when 
P = or 1. Moreover, based on more physical arguments and following Kronmiiller's analysis,— 
Drobrynin et al.,— derived an expression for the effective demagnetizing field which contains a 
term that depends on the difference of the outer and inner demagnetizing factors as in Eq. (fTTT ). 
which provides support for Netzelmann's interpolation.— 
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The expression obtained for the effective demagnetizing field in the saturated state, Eq (UOl ). 
shows that this differs from that of the homogeneous material defined by the outer volume when 
P < 1. For the particular case of a perpendicularly magnetized assembly with an outer volume of 
a thin film, as shown by Eq. (I22l) . the effective demagnetizing factor is lower than the continuous 
film, < in provided that P < 1.^ 

The extensions to include the configuration dependence are done by taking the product between 
the interaction term for the saturated case and Eq. (|T6l l on the component of the demagnetizing 
field along the easy axis, Eq. (fTSi) . and on the total anisotropy field, Eq. (|2T1) . 

An important result is that the interaction field can be expressed different forms each with 
different physical meaning, which however, are related among them. The dipolar interaction con- 
tributes to the components of the demagnetizing field, Eq (flOl) . as well as to the total magnetostatic 
anisotropy field, Eq. (fT2l) . Then each of these contributions can be extended to non saturated states 
assuming an interaction field of the form am, where the coefficient a is given by Eq. (fT9l) . The 
FMR and magnetometry measurements done on arrays of nanowires provide a validation of this 
approach, while showing that different measuring techniques can provide values that correspond 
to different forms of expressing of the interaction field. 

From the extension done to include the configuration dependency on the interaction field, the 
hysteresis loop shearing can be attributed to the magnetization dependent component of the effec- 



tive demagnetizing field along the easy axis. As discussed in Ref. |23|, assuming an interaction field 
of the form am, the shearing can be corrected once the value of the coefficient is known, which; 
as shown here, is given by Eq. (fT9] ). 

The coercive field in an assembly is defined with respect to the total or effective anisotropy of 
the entire system, regardless of the specific internal composition, and in the Stoner - Wohlfarth 
model, these quantities are equal. If we call Hc{0) = MgAN the coercive field of the isolated 
particle (P = 0), rearranging Eq. (fT2] |. the coercive field of the assembly is given by, 

H, = H,{0){1 -P) + MsAN+P, (36) 

which depends on the packing fraction as well as on the outer volume. The first term is the well 
known coercive field dependence on the packing fraction in assemblies where the shape anisotropy 
prevails The second term represents the outer volume anisotropy contribution, which implies 
that the coercivity of the assembly is shape dependent. So for a given assembly with fixed packing 
fraction, the coercive field will change between samples with different outer volume, as recently 
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reported for granular hard magnetic films. — Moreover, for the particular case when the outer vol- 
ume is isotropic and AA^+ = 0, the second term vanishes. To further emphasize the importance of 
the second term in Eq. (|36] ). notice that if only the first term is considered, the coercivity always 
decreases when the packing fraction increases, regardless of how the particles are distributed .-^241 
This behavior is only consistent with the case where the interaction between the particles is anti- 
ferromagnetic. Indeed, for example, consider the linear chain of cylinders discussed above, each 
with their easy axis parallel to the chain axis, the interaction is ferromagnetic and as the packing 
fraction increases, the coercivity of the entire chain should increase from the value of the isolated 
cylinder of finite aspect ratio to that of the infinite cylinder. This behavior is described by Eq. (|36] | 
when the second term is taken into account. 

Since the total anisotropy field, Eq. (fT2] ). depends explicitly on this outer shape anisotropy 
MsAN~^P, then variations on the shape of the outer volume results in changes in the demagnetiz- 
ing field, the anisotropy and, as just mentioned, in the coercivity. This contribution is associated 
only to the dipolar interaction, so at low packing it is weak, while at higher packing fractions it 
can induce appreciable effects, and can be used as an additional parameter to control the magnetic 
properties of the assembly. An interesting case is when the assembly is made by spherical parti- 
cles, since their shape anisotropy vanishes and the effective anisotropy field is given only by the 
outer shape anisotropy. Then a chain of spheres will behave close to a cylinder (its outer volume) 
and the entire chain has a uniaxial anisotropy whose origin is the interparticle interaction field,— 
while the same particles forming an spherical assembly will be isotropic. Although this is a well 
known result, it shows that this effect follows from the key role played by the outer demagnetizing 
factor on the value of the interaction field, which recently has been related to novel effects in gran- 
ular hard magnets,'^ composites,— and dipolar ferromagnetic order in superparamagnetic particle 
monolayers.— 

Another effect related to this contribution is the magnetization reorientation transition induced 
by the interaction field, which was discussed for the linear chain of cylinders, but which has also 
been reported for other assemblies.— -i^i^^ As discussed above, this transition results with in- 
creasing the packing fraction when the shape anisotropy of the single particle and the dipolar term 
have opposite signs. Qualitatively, this can also be interpreted as a competition between the shape 
anisotropics of the inner and outer volume and how their respective easy axis are oriented with 
respect to each other. At very low packing, the easy axis of a given particle in the assembly is that 
of the inner volume, while as P — )■ 1 the easy axis corresponds to that of the outer volume. If 
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both inner and outer volume have their easy axis in the same direction, there will be no reorienta- 
tion transition of the magnetization, but if their easy axes are not aligned, this transition will take 
place above a certain critical packing fraction. In this sense, both inner and outer volumes should 
be anisotropic, since if either one of them is isotropic, there will be no transition. Moreover, an 
anisotropic outer volume is a necessary condition to induce this transition as deduced from Eq. 
(l36l) . Noting that in this expression He is the total shape anisotropy and Hc{0) = MgAN corre- 
sponds to the inner volume shape anisotropy, it follows that the total anisotropy will never change 
sign if the second term is not included. 

The validity of this mean field approach is determined mainly by the extent to which the indi- 
vidual particles in the assembly remain bistable and are able to rotate independently from other 
particles, this is, that no collective effects take place while individual reverse their magnetic state. 
Starting with a very diluted assembly of bistable particles, both conditions are fulfilled, however 
as the packing fraction increases, each particle feels the interaction field which increases as the 
distance between particles is reduced. In this sense both bi-stability and the ability to reverse inde- 
pendently are susceptible to change as a consequence of any finite value of the interaction field. On 
one hand, bi-stability can be related to the strength of the coercive field of the individual particle, 
or to the height of the energy barrier that separates the two minima configurations (+m and — m). 
On the other hand, the ability of the particles to rotate independently also depends on the value of 
the coercive field and the width of the intrinsic switching field distribution. So in both cases, the 
competition between the height of the intrinsic energy barrier of each particle and the magnitude 
of the interaction field determines if these conditions are fulfilled and in this sense, there is no 
unique criteria but rather it will depend on the type of material and the main contributions to the 
energy. 

In particular, the height of the energy barrier of the particles depends on the magnetocrys- 
talline anisotropy. Since this contribution is intrinsic to each particle, it follows that the effects of 
the dipolar interaction are expected to become less relevant as the magnetocrystalline anisotropy 
of the individual particles increases.—^ In this sense, the magnetic hardness parameter— k = 
[li^Mcl/lSvrM^)]^/^, provides a useful empirical parameter of the height of the energy barrier. 
In particular, the higher the value of n, the higher the packing fraction attainable without loss of 
bistability or independent reversal. On the other limit, this is, low values of k both the bi stability 
and independent reversal are expected to be limited to low and moderate packing fractions, as 
suggested by previous reports which show that magnetic percolation in assemblies takes place for 
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packing fractions between 20 and 40%.-^2i^ 

V. CONCLUSION 

In conclusion, mean field expressions have been derived for the demagnetizing field and the 
magnetostatic anisotropy field which include the magnetization dependent part of the interaction 
field for assemblies of magnetic particles with arbitrary inner and outer demagnetizing factors. 
Special emphasis was given to 2D arrays of nanomagnets with perpendicular magnetization, where 
the expressions for the demagnetizing field and total anisotropy field have been derived. The 
model was successfully tested using FMR and magnetometry measurements done on arrays of 
cylindrical nanowires, where the equivalence between the interaction field determined by each of 
these techniques was shown. This formalism provides a simple mean field framework to describe 
magnetostatic effects in a wide variety of magnetic assemblies and composites. Moreover, it shows 
the role played by the outer demagnetizing factor in the value and characteristics of the interaction 
field which could be used to fine tune the magnetic properties of the assembly or composite. 
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